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A key lesson of the decoherence program is that information flowing out from an open system 
is stored in the quantum state of the surroundings. Simultaneously, quantum measurement theory 
shows that the evolution of any open system when its environment is measured is nonlinear and leads 
to pure states conditioned on the measurement record. Here we report the discovery of a fundamental 
relation between measurement and entanglement which is characteristic of this scenario. It takes 
the form of a scaling law between the amount of entanglement in the conditional state of the 
system and the probabilities of the experimental outcomes obtained from measuring the state of the 
environment. Using the scaling, we construct the distribution of entanglement over the ensemble 
of experimental outcomes for standard models with one open channel and provide rigorous results 
on finite-time disentanglement in systems coupled to non-Markovian baths. The scaling allows the 
direct experimental detection and quantification of entanglement in conditional states of a large 
class of open systems by quantum tomography of the bath. 

PACS numbers: 03.65.Ta, 03.65.Ud, 03.65.Yz, 03.67.Bg, 03,67,Pp 



Arguable, a denning feature of a real quantum system 
is the nature of its unavoidable interactions with its sur- 
roundings. It is well known that they are determinant of 
the manners the system can evolve, being even respon- 
sible, in occasions, of the total loss of the system quan- 
tum character pQ. Less discussed, however, is the fact 
that they also dictate the ways a system can be observed 
and information can be extracted from it. Thus, for in- 
stance, a system can be indirectly probed by perform- 
ing a projective measurement on its environment. Ignore 
the outcome of this measurement, and the information 
lost implies that a pure state of the system prior ob- 
servation is transformed into a statistical mixture [5J [3] . 
If on the other hand, the result of the measurement is 
physically stored with full fidelity the state of the cen- 
tral system immediately after measurement is again pure, 
yet conditioned on the particular, stochastic outcome: 
the state, representing our knowledge about the system, 
"collapses" [H [5]. Conditional states not only have an 
intrinsic importance as the outcome of physical measure- 
ments. They provide also an exact representation of the 
unconditional dynamics, with implications for the field 
of quantum information [SHE] where characterizing the 
amount and kind of non-classical correlations in uncon- 
ditional states of open systems remains a main theoretical 
challenge [Hi [10] . In this order of ideas, a proper under- 
standing of the interplay between measurements on the 
bath and entanglement properties of conditional states, 
to our knowledge a problem not yet explored, is clearly 
necessary. 

Under the only assumption of a strictly local dynamics, 
in the present article we derive a scaling law relating the 
probability distribution of experimental outcomes from 
measurements on the environment, as given by the Born 



rule, to the amount of multipartite entanglement of the 
pure conditional state of an open quantum system. As we 
demonstrate, our findings are universal: They are inde- 
pendent of the number and dimensions of the constituent 
subsystems and hold for any of the so-called G-invariant 
measures of entanglement dHHl], e.g., Wootter's concur- 
rence for two-qubit systems [T^, characterized by their 
invariance under local invertible transformations. No- 
tably, these results allow for a characterization of multi- 
partite entanglement in the presence of both Markovian 
and non-Markovian environments, thus making contact 
with state-of-the-art experimental techniques |14j . Fur- 
thermore, we show that our scaling law suggests a proto- 
col for a direct measurement of entanglement in a large 
class of systems, including all Markovian and all non- 
demolition baths: As it turn out, the inverse numerical 
value of the Husimi function of the bath is a measure of 
the amount of entanglement in the conditional state. 

Our work is partially motivated by recent results in 
the study of non-classical correlations in unconditional 
states, where a conservative approach based on explicit 
calculation of the entanglement measures in the uncondi- 
tional mixed state faces tremendous difficulties due to an 
intrinsic high dimensional optimization problem |15[ 116] . 
The severity of this hurdle has motivated intense re- 
search in the general (algebraic) aspects of the interplay 
between unconditional time evolution and entanglement 
measures. An important step in this direction was the 
discovery of a scaling property of the G-invariant mea- 
sures |1U II 7] , and its experimental verification for two- 
qubit systems with one open channel [18] , 

However, the scaling property of unconditional states 
critically depends on the linearity of the evolution [11] 
through its Kraus representation [9 . Therefore, an ex- 
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tension of such algebraic properties to the conditional 
case with its intrinsic nonlinear dependence on the initial 
state requires both different mathematical methods and 
physical concepts. Our results fill this gap at a moment 
when there is an increasing interest on schemes with a 
less restrictive role for the environment, where the flow 
of information from the system to the bath takes an ac- 
tive part [6HH1 H9H23] . 

Open systems in photonic environments 

We consider a central-system-plus-bath approach [U [5] 
where a multipartite central system (CS) made of N sub- 
systems with dimensions dS 1 ' , . . . , dS N ) is coupled to a 
collection of independent photonic baths (channels) act- 
ing locally on each subsystem. If initially pure, after 
transformation to the interaction picture, the state of 
the combined system is given at time t by 

|*(t)) = Te~^ ti {6as{s)+6l(s))ds \<t>{$)) ® |0). (1) 

In the above equation |</>(0)) is the initial state of the CS 
and |0) the ground state of the bath, T stands for the 
time ordering operator, Hcs(s) is the CS Hamiltonian 
acting locally on the subsystems, and 

fe : A 

accounts for the interaction between CS and bath 
through local traceless operators with gk,x the cou- 
pling strength of the system and the A-th mode to the 
fc-th channel, where photons of frequency u>k,x are created 
or annihilated by a\ x and (ik,x, respectively. Memory ef- 
fects in this approach are characterized by the correlation 
functions 

/>oo 

«*(*) = E \9k,x\ 2e ~ iUk ' xS = / h(u)e- ius du, 
x Ja 

where Ik{oS) is the spectral density of the bath's fc-th 
channel (see Supplementary Information). 

Universal scaling of multipartite entanglement 

Following [24H26] . the total state in equation 

can be always written as 

|*(t)) = Z- 1 J e-^|V(a, t)) ® \a)da, (2) 

where \tf>(a,t)) is the (unnormalized) state of the system 
relative to the bath (normalized) coherent state \a), Z 
is a normalization constant, which will drop out from all 
the final results, and we define da : = ]X > dtftak^xd^ak x 
with \a\ 2 = J2k,x \ a k,x\ 2 - 



If at time t > a generalized measurement of the bath 
is performed, the probability P(a,t) that we obtain \a) 
is [31 [25] 

= WM)|#M)>, (3) 

with P(a, 0) = Z~ 1 e~^ 2 , and the state of the central 
system collapses to the conditional state 

Ma,t)> = WM)^M)}- 1/2 |^M)) ■ (4) 

To proceed further, we notice that the unnormalized 
conditional state \4>{a,t)) can be related to the physical 
initial state of the CS by means of [24] 

\i)(a,t))=U(a,t)\<t>(0)) (5) 

with 

U(a,t) = Teti[-% 6 ™ + ^k( z ^ J *- J t 6 <>( s d} ds , 

a non-unitary propagator depending on the measurement 
outcomes {a} through Zk(s) = J2x 9k,xa-k,xe~ lUk,xS and 
Ok{s) = Jq a k (s - s')Ok{ak,s,s')ds'. Although in gen- 
eral the operators 0&(cife, s, s') are defined only pertur- 
batively for each particular choice of ife(w) and exact 
closed expressions are available for several cases of phys- 
ical interest [24H26] . 

Having a linear relation between \t()(a,i)) and |</>(0)), 
we can turn to the main issue of this work: the amount 
of entanglement G(\<p(a,t))) in the conditional state 
\4>{a,t)) as quantified by a G- invariant measure. As we 
now show, the answer to this question follows smoothly 
from two defining properties of G (i) It is invari- 

ant under local linear transformations with unit determi- 
nant, (ii) It is homogeneous of degree two, i.e., G(u\(j>)) = 
|u| 2 G(|</>)) for all complex u. Indeed, from property (ii) 
and Eqs. (|4|,([5]), the entanglement of \<f>(a,t)) can be 
written as 

G(\<j>(a,t))) = MM)Ma,0> _1 G(«M)l#>)>)- 

We now use property (i), in conjunction with equa- 
tion ([3|, to obtain the central result of this paper, 

where the scaling function 

f(a,t) = e~^ k <Xk(s-s')Tr[jld k (a k ,s,s')]dsds' 

is independent of the functional form of G(\4>)); a con- 
sequence of the use of the identity detTe^o A ( s ) ds = 
e/o T ? A ( s ) ds anc j the fact that the Hamiltonians and baths 
act locally in the subsystems (see Methods section). 
Equation ([6]) provides us with an explicit scaling between 
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the amount of entanglement in the conditional state of 
the central system and the probability distribution of the 
outcomes of the measurement on the bath. It expresses 
universality of the multipartite entanglement in the sense 
that for a given outcome of the measurement, all normal- 
ized G-invariant measures of multipartite entanglement 
in the conditional state have exactly the same numerical 
value. Within the standard system-plus-bath approach 
for local systems, our equation ([6]) is completely general. 

To illustrate the relevance of scaling law ^ , in the re- 
mainder of the paper we use it to address some crucial 
questions in the actual theory of multipartite entangle- 
ment. 



Distribution of multipartite entanglement 

Perhaps the most pertinent question at this point is: 
what is the probability of having a given amount of entan- 
glement in a system just after its photonic environment 
has been measured? Let us remark that the common 
approach of using uniform distributions for \(f>) |27[ 128] 
does not hold here. Once the measurement process is 
specified, i.e., a choice for the measurement outcomes 
is made, unitary invariance as a guiding principle is re- 
placed by the specific quantum-mechanical distribution 
of outcomes P(a,t), as given by the Born rule. Accord- 
ingly, the distribution Pc(x,t) of the normalized entan- 
glement x(a, t) in the physical ensemble of pure and ran- 
dom conditional states of the CS must be calculated from 
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FIG. 1: Probability distribution of entanglement. 

Probability distribution of scaled multipartite entanglement 
in the CS over the ensemble of experimental outcomes {a} 
of an amplitude damping bath for different initial states: a, 
pn = 1/2, pio = 0, b, pn = 3/4, pio = 0, c, pn = 1/2, 
pio = 1/4, d, pn = 1/2, pio = 1/3. The figure is the same for 
any G-invariant measure of a system of qudits if only a single 
qubit without internal dynamics is coupled to a zero temper- 
ature Markovian bath. The black continuous line shows the 
maximum possible value of entanglement a; max (p) while the 
red dashed line corresponds to the mean entanglement x(p). 



Pc(x,t) = P(a,t)S (x — x(a,t)) da . 



(7) 



which can be evaluated using the scaling rule ^ and is 
the same for all G-invariant measures. In particular, its 
first moment, i.e., the mean entanglement x(t), is given 
by 



P G (x,t)xdx= P(a,0)f(a,t)da. (8) 



x(t) 



To exemplify our argument we consider a specific and 
experimentally relevant case |18) : a multipartite CS 
where only one qubit is coupled to a single open chan- 
nel (k ="ch"). Hence, the information about the initial 
state |</>(0)) is fully encoded in the initial reduced density 
matrix of the coupled qubit, 



Pch(0) = 



Pn Pio 
Pio Poo 



Let us start by considering a Markovian bath at zero 
temperature where J = ■Jjfr- = 0(a,s,s') 25 , with 
<7_ the two-level deexcitation operator, 7 the decay rate, 
and I(u>) = 7T _1 . In order to focus on decoherence ef- 
fects, we assume frozen internal dynamics. In Fig. [T] 
we show Pc(x,p) in terms of p c h(0) and the scaled time 



P(t) 



-it 



(see Supplementary Information). In all 



cases we observe that the entanglement distribution be- 
comes broader as time increases to finally collapse asymp- 
totically when entanglement from CS disappears. The 
sharp boundary exhibit by Pc(x,p) (black continuous 
line) is given by the condition 



V 1 -p f 1 - pnP 



Pll 



P11 
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Thus, for maximally entangled states, i.e., states with 
Tr/5^ h (0) = 1/2, the maximum possible value of entan- 
glement x max (p) is strictly decreasing with x max (p) < 1, 
as can be seen in panel la. However, for initial states with 
T r Pch(0) > 1/2; there are regions in which x max (p) can 
increase even up to a; max (p) > 1, as displayed in panel lb. 
This behaviour should be contrasted with the evolution 



of the averaged normalized entanglement x(p) = y/l — p, 
which is independent of the initial state and strictly de- 
creasing (red dashed line). 

A quite different picture arises when the channel cor- 
responds to a non-demolition bath. Thus H\ 



(ch) 
CS 



oc J = 



^a z = 0(a, s, s 1 ) [25] , tr z is the third of the Pauli's ma- 
trices, and A 2 the "phase-flip" rate. In this so-called de- 
phasing channel, the distribution Pg{x, t) is independent 
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FIG. 2: Probability distribution of entanglement. 

Probability distribution of scaled multipartite entanglement 
in the CS over the ensemble of experimental outcomes {a} of a 
non-demolition bath for different initial states: a, pu — 1/2, 
b, pn = 1/4. The figure is the same for any G-invariant mea- 
sure of a system of qudits if only a single qubit is coupled to 
a non-demolition bath. The black continuous line shows the 
maximum possible value of entanglement a; max (p) while the 
red dashed line corresponds to the mean entanglement x(p). 
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FIG. 3: Evolution of entanglement distribution and 
mean entanglement, a, Probability distribution of multi- 
partite entanglement at different values of the scaled time p. 
The black dots mark the mean entanglement x(p). b, Loga- 
rithm of the average entanglement X £LS cL function of A 2 t for 
ohmic (continuous lines) and superohmic (dashed lines) envi- 
ronments and different values of the effective Debye frequency 
(lu = cjd/A 2 ). Both behaviors demonstrate a transition from 
exponential to much slower algebraic decay, also characteristic 
of the approach p(t — > oo) — > 1. 



of pio and has a sharp, time independent boundary given 
by x max = (4pnpoo) • For maximally mixed p c h(0) 
the upper limit is £ max = 1, while x max has square-root 
divergences for pure initial states pu = 1,0. In Fig. [2] 
we plot Pc(x,p) for the Markovian case with = n^ 1 
and rescaled time pit) = 1 — c~ A *. 

For this class of environments the maximum of 
Pg{x, t), which signals the most likely values of the mul- 
tipartite entanglement, do not lie near the average value 
x(t). Instead, the distributions peaks around the bound- 
ary x = x max , where it has an integrable divergence (See 
Fig. |2]). This is clearly seen in the Markovian case with 
Pu = 1/2, where 



Pc(x,p) 



2(1 -p) 



1 



7r|log(l -p) \ X \/l - x 2 
log 2 ( l+VIEl 



x exp 



2|log(l-p)| 



and that we show in Fig. [3j When time increases a sec- 
ond maximum appears near (but not exactly at) x = 0, 
showing that the average x(p) (red dashed line in Fig. [2] 
and points in Fig. [3]) misses the most likely values of 
multipartite entanglement. 

In this case, the dependence of our results on the mem- 
ory scale of the environment is easily checked. For this 
we introduce non-Markovian effects using a spectral den- 
sity I(w) 7^ 7T _1 [3] (see Supplementary Information), 
and find that this amounts only to a redefinition of time 
t — > q(t) in the purely ohmic result by means of 



?(*) = 



J( W ) 



sin ujt/2 
cot/2 



did . 



The parameterization p(t) = 1 — e^ 2 ^*) for the time 
is of most convenience in this case, leading to the mean 
entanglement x(p) = ^/l — p. The function p(t) rules 
the actual speed at which the distribution and its first 
moment attain their asymptotic form. It is given for 
exponential cut-off at Debye frequency uid by 



p(t » io d x ) = 



(l+^2)l/™ £ 

l/(l+o;2t 2 ) 1 /™ 



-A 2 t 



(ohmic) , 
(superohmic) , 



where we introduced u) = cj^/A 2 . The above relation im- 
plies that the decay of the averaged entanglement is only 
algebraic for superohmic environments (see Fig. [3|. Re- 
markably, however, the approach to the limit u^t — > oo 
is such that for dephasing channels the average entangle- 
ment is never zero for finite times. 



Finite disentanglement time 

As discussed above, the most likely values of mutipar- 
tite entanglement in conditional states are in general not 
related to their average value. A prime exception takes 
place when there exists a finite time t = t such that 
x(t) = 0, since it entails that at r all entanglement mea- 
sures of the multipartite conditional states vanish identi- 
cally over the whole ensemble of experimental outcomes. 
The conditions for this event to occur as well as the de- 
pendence of the disentanglement time t with the parame- 
ters of the system are issues of obvious practical interest. 
Here we answer these questions for a large class of sys- 
tems where at least one channel k = I satisfies 



Oi(ai,s, a') = Oiui(s, s') 



(9) 
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for some function ui(s,s r ), and operator O/. This class 
of systems includes all Markovian and all non-demolition 
baths, the Jaynes-Cummings model and the damped har- 
monic oscillator among others |25j . 

We notice two immediate consequences of the form of 
equation ^ for the mean entanglement: First, if |9| 
holds for all channels, the mean entanglement reduces to 



x(t) = c 



_ Q "E* ^Tr\jl6 h ]fZ f° a h ( s - S ')u k ( s , S ')dsds' 



(10) 



Second, in the Markovian case, u(s,s') = 1, O = J, 
x(t) is independent of the internal dynamics and strictly 
exponential; no finite r exists and the mean entanglement 
disappears only asymptotically in time. 

More generally, the product form of both P(a, 0) and 
/(a, t) implies that the vanishing of the factor associated 
with a single channel will automatically disentangle the 
whole system. Let us focus on this particular channel 
and assume that it couples a qubit with frozen inter- 
nal dynamics to a non-Markovian zero-temperature bath 
with memory kernel a(s) — (ujd/2)e~ u ' d \ s \ . The exact 
dynamics of the total system for arbitrary a(s) is ob- 
tained with the choice O — ^/yct- as long as the func- 
tion u(s, s') satisfies u(s,s') = c(s)/c(s') and dc(s)/ds = 
—7 J* a(s — s')c(s')ds' , which can be solved for the case 
we are considering (see Supplementary Information) . Af- 



ter substituting our result into ( 10 ) we obtain 



x(i) oc exp < —7 



|g sinh fls 



cosh ils 



ds 



where fl — (wd/2) \/l — p and \x = 2j/wd- We notice 
that for weak and strong coupling (/1 < 1) the average 
entanglement is not zero for all finite times. In contrast, 
in the ultrastrong coupling regime (/1 > 1) we get x(t) = 
for 



7T 



arctan \J p — 1 



rigorously showing the existence of finite disentanglement 
times in all conditional states of arbitrary multipartite 
systems where at least one single qubit is coupled with 
a non-Markovian bath. Remarkably, the disentangle- 
ment time increases with p, for strong enough coupling, 
as shown in Fig. |4j 

Our last considerations are directly extended to the un- 
conditional dynamics of the open system (entanglement 
sudden death [29 -34 ). The mixed state describing the 
unconditional evolution is given in complete generality 
by an average of the pure conditional states [2~4T - f2l)] 

P unc (t) = J P{a,t)\<f>{a t t)){4>{a,t)\da. 

Therefore, the set {|0(a, i)), P(a, t)da} constitutes a valid 
convex decomposition of p unc (t), and since the amount 




FIG. 4: Disentanglement times, a, Dependence of the dis- 
entanglement times r with the effective coupling n — 27 /w^ 
for an arbitrary multipartite system where a qubit without 
internal dynamics is coupled to a zero-temperature bath with 
a characteristic memory scale u7 . b, Same as in a showing 
that for large coupling the natural scale for the disentangle- 
ment time r is given by (-yuid) -1 ^ 2 ■ 



of entanglement G(p unc (t)) is defined by a convex roof 
construction, we have 



G(p u 
G(p» 



= (0)) 



< x(t). 



Therefore all measures of multipartite entanglement of 
an open system interacting with an arbitrary photonic 
environment must be zero at a finite time if at least one 
qubit is strongly coupled with a non-Markovian bath at 
zero temperature. The important extension of this result 
to finite-temperature baths is presently under study [35] . 



Tomography of entanglement 

To conclude, in this last section we show how scaling 
law ^ provides us with a geometrical and operational 
interpretation of the amount of entanglement in condi- 
tional states of multipartite systems. This follows from 
the observation that, associated with the pure total state 
IvE'(i)), the unconditional density matrix describing the 
environment p E (t) = Tr CS | ^(t))(^f(t) \ [3S] has the fol- 
lowing representation in coherent states 

Q(a,t)=Z- 1 (a\ p&cs|*(t)><*(t)|]|a>, 

where Q(a, t) is the so-called Husimi function of the bath 
[37] . Using equation |2]) then leads to 

Q{a,t) = P(a,t), 

and therefore equation (|6| relates the possible values of 
the entanglement of the conditional state with the mor- 
phology of the Husimi function of the bath; a quantity 
accessible to direct measurement. 
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In the following we will consider systems where equa- 
tion ([9]) holds. In this case we easily obtain 



x{a,t) = _ (Na _ Nnl) Q{a',t) 
x{a',t) 6 Q(M)' 



(11) 



where N a = \a\ 2 is the mean number of photons in the 
state \a). Hence, the numerical value of the entanglement 
of the conditional state after a single measurement with 
experimental outcome {a'} together with the uncondi- 
tional Husimi function of the bath Q(a' , t) automatically 
predicts all the entanglement measures of all possible 
multipartite conditional states. Quantum tomography of 
the bath is actually measuring entanglement of the con- 
ditional states. 

Finally, known upper bounds for the Husimi function 
[57] can be used to provide a lower bound to entangle- 
ment 

x(a,t) > e~ Na x(t), 

which is particularly useful in the experimentally acces- 
sible case where \a) lies close to the ground state. 



Summary 

The derivation of a scaling law relating the amount of 
entanglement in the state of a multipartite system just 
after its photonic environment is measured, to the prob- 
ability distribution of the experimental outcomes, consti- 
tutes the central result of this paper. It indicates that 
for bosonic environments, multipartite entanglement is 
severely constrained by indirect measurement. Specifi- 
cally, we use equation ^ to demonstrate that all pos- 
sible normalized measures of multipartite entanglement 
in the conditional state have the same distribution over 
the ensemble of measurement outcomes of the bath. We 
construct the full distribution in standard cases with one 
open channel and show that for non-demolition environ- 
ments non-Markovian effects appear only through a re- 
definition of the time variable. We rigorously prove dis- 
entanglement in both pure conditional and mixed uncon- 
ditional states of multipartite systems where one qubit is 
coupled with a zero temperature non-Markovian bath, 
with a counter intuitive dependence of the disentangle- 
ment times on the coupling strength. Finally, we showed 
that due to the scaling, for a large class of systems, multi- 
partite entanglement of the conditional states is directly 
measured by quantum tomography of the bath through 
the Husimi function. 



Methods 

Scaling function. Here we calculate explicitly the 
scaling function in equation (JsJ) . By definition, f{a,t) 



must satisfy 

G(U(a,t)\<j>(0))) = f(a,t)G(\m)) 

with 

U(a,t) = Te^1-^ cs+ ^( 2 * (s)ife ~ J fc 6fc(s) )K 

Assuming that the Hamiltonian -ffcSj the Linblad op- 
erators Jfe, and the convolutionless operators Ok{s) act 
locally on the subsystems, the operator U admits the fac- 
torization 

N 



W(a,f) = []wW(fl«,(). 



Here 



W ( a W , t) = Te f ° (<m -^ hW & h w ( s ))] <■ 



U 



is a local operators acting on the ith subsystem, and 
{fc^} denotes the indexes of the Lindblad operators act- 
ing on the same subsystem. We now introduce operators 

U^(a^,t) = |det wWtaW.tJI-^WWfcW,*) 

with |det 17 W| = 1, since is the dimension of the 
Hilbert space of the ith subsystem, and use the homo- 
geneity of the function G to obtain 



N 



2/d 1 -' 



G(U{a, t)\(j>{0))) = Y[ det U {t) (a« , t) 

i=l 

/ r N 

IJl7« (a«,t) 



x G 



i=l 



10(0)) 



Since by definition G(\4>}) is invariant under local trans- 
formations with determinant one, we can identify 



N 



f(a,t) = Y[ det W (i) (fl (i) ,t) 



2/S 1 



To proceed we focus on the individual terms in the prod- 
uct. Under mild conditions, the determinant of U^ 1 ' ad- 
mits the representation 

det W«(a«,,s)=e Trl °s * t0 <° (< H 



and therefore 



log 



det W (l) (a W ,s) 



2/d 1 ' 



23* 



Tr logtfW(aW,«). 



To evaluate this last expression we start from the equa- 
tion of motion for WW, 



9 m(a®,s) = 



ds 



Lfc(« 



U {i) {a (i) , 
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to write 

d 



ds 

5tTr 



log 



det U w {a^,a) 



2/d« 



^s + E(^coW^w-^(oO fc «W) 



fc(<) 



where we used the cyclic property of the trace. Further 
simplifications take place after noticing that Hcs is her- 
mitian and that by construction </& is traceless, then 



d_ 



log 



det U (i) {a^,s) 



Tr ^Jt (i) fc(i) ( s ), 



which is trivially integrated. Substituting the result 
in f(a,t) and collecting the sums in the exponents via 

E 4 Efcw = Efe wc nnall y g et 

/(a,t)^c-^^^ Tr [ J ^ (s) ] ds , 
thus completing the construction of the scaling function. 
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SUPPLEMENTARY INFORMATION 
Open system Hamiltonian 

In this section we derive within the usual system-plus- 
bath approach (see Refs. 4, 24, 25 and 26) the Hamil- 
tonian of an open system appearing in equation (1) in 
the main text. For this, we consider the system dynam- 
ics in a convenient interaction picture and show how to 
redefine the Lindblad operators to make them traceless. 
For simplicity we consider a single channel for the deriva- 
tion; the generalization to arbitrary number of channels 
is straightforward and will be presented at the end of the 
section. 

Our starting point is the total Hamiltonian H describ- 
ing a central system (CS) linearly coupled with a bosonic 
bath (B), 



H = H 



cs 



Hi 



CSB 



Hp 



(12) 



with environment and coupling terms given by 
H B = ^ w\a\a x , 



H 



CSB 



Lot - L ] 



We consider a bosonic bath with frequency uj\ in its Ath 
mode, for which a\ and a\ correspond to anihilation and 
creation operators, respectively. The operator L acts on 
the CS and is usually dubbed Lindblad operator. To- 
gether with the real coupling constants g\ it describes the 
interaction of the central system with all the modes of the 
bath. The total Hilbert space H = Ties ® Hb of the sys- 
tem is spanned by product states | ) <g> \ai, . . . , ax, ■ • •), 
where {|ej}} is any basis of the CS Hilbert space, and 
\a) = \a%, . . . , a\, . . .) is a normalized coherent state of the 
environment, i.e., a common eigenstate of the annihila- 
tion operators of the bath, d\\a) = a\\a) with (a\a) = 1. 
Introduce now the transformations 

t .9\ 

a x = b\- i — a , 

L = J + a, 

with complex constant a. Observe that the coherent 
states | a) are also eigenstates of b. Then, up to a physi- 



cally irrelevant constant term, Hamiltonian ( 12 ) becomes 



H — Has — 



a*J + ajA + 



E 

A 



w\b\b 



A OA 



In order to obtain a traceless coupling operator J, we take 
advantage of the freedom introduced by a and choose 
a = dr lr Tr L, where d is the dimension of the Hilbert 
space of the CS. Thus the hermitian redefinition of the 



CS Hamiltonian 

2 

H cs ^H C s-J2^{ a *J + a ^) ( 13 ) 



allow us to work with traceless Lindblad operators in 
full generality. Since this change, up to a constant, is 
the only consequence of the transformation, we stick to 
the initial notation and write the total Hamiltonian as in 



equation (12 1 but with J instead of L, and keep in mind 



renormalization ( 13 1 for the CS Hamiltonian. 



The interaction picture for the total time-evolution op- 
erator U(t) = e~i Ht is defined via the unitary transfor- 
mation 



#int(*) =' 



U(t)e 



and absorbs the explicit dependence on the Hamilto- 
nian of the bath into a redefinition of the bath operators 
^a^W — a\e~ lWxt , rendering the transformed Hamilto- 
nian explicitly time-dependent. Formal solution of the 
Schrodinger equation in the interaction picture yields 

where Hcs( s ) is the renormalized CS Hamiltonian, T is 
the time-ordering operator, and 



Hi(s) = i^5fe,A (jkd{ 



fc,A 



is the interaction Hamiltonian. In this last step we al- 
ready introduced the result for the general case with an 
arbitrary number of channels indexed by k. 



Entanglement distributions 

Here we evaluate the probability distributions of mul- 
tipartite entanglement for conditional states in the CS 
over the ensemble of experimental outcomes of measure- 
ments on the environments studied in the main text. Us- 
ing the scale law stated in equation (6) of the main text, 
the probability distribution of entanglement, as given by 
equation (7) of the main text, takes the manifestly uni- 
versal form 

2 J x 



xS x „ )da, 

\ Tr p CS (Op{a,t)iU(a,t) J 

where we used equations (3) and (5) from the main text 
for simplifications. The particular dependence of the 
scaling function f(a,t) and U(a,tyU(a,t) on the bath 
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measurement outcomes {a} must be found in each par- 
ticular case of interest from their definitions (see Methods 
section in the main text): 

/(a t) = C~^ fe & fa a k (s-s')Tr[Jld k (a k ,s,s')}dsds' 

U(a,t) = Teti[-^ Aca+ ^( z ^ s)Jk - J i 6k ^)] ds . 

Further simplifications in the expression for the en- 
tanglement distribution are obtained if we assume that 
only a single subsystem is coupled to a channel and 
use the locality of Hcs- In this case, the part of 
Tr pes (0)W (a, t)^U (a, t) which traces over the subsystems 
not coupled with the open channel can be calculated, 
leaving only the reduced density matrix describing the 
initial state of the coupled subsystem (denoted by "ch"). 
If in addition we suppose that 

for some operator O and function u(s, s'), then /(a, t) — 
fit) is independent of {a}, and all variables a,(k^ch) in 
the probability distribution can be integrated out by per- 
forming the Gaussian integral left for them. All together, 
the key technical difficulties we must face to obtain the 
entanglement distribution are the evaluation of the func- 
tion 

F(a (ch) ,t) = Tr ch p ch (0)ti^ (a (ch) , i)tw( ch >(a (ch) , i), 
and the calculation of the integral 



P G (x,t) 



fit) 



3 -k(ch)l 



-S x 



-(ch) 



/(*) , , 

(14) 

We now follow this program for the two physically rele- 
vant situations with one open channel considered in the 
main text. 



the bath defining the open channel. The above equation 
can be solved by introducing the interaction picture in 
ri_ and exploiting the commutator [ri-,(T±] = ±<7 ± . the 
result is 

U(a,t) = e -i^-t e V^fo z '( a ^ i3d ^-. 

Explicit calculation using the algebraic properties of the 
Pauli matrices then gives 

U(a t tyti(a,t) = 7 + (e" 7 * - 1 + l^jftt)! 2 ) n_ 
+ [a T g{t)°+ + (cTftt))'*-] , 

where 



a r 9{t) = \/ll z*(a,s)e * s ds 
Jo 

e -i(«A-ii)t _ ! 

A 



Consequently, F(a, t) depends on the integration vari- 
ables a\ only through the linear combination y(a,t) — 
a T g(t). The distribution function for y — y(a,t) is read- 
ily calculated, 



./-',(.//./)= / ^-5 [>, -- u T <][!)) d„ 



1 



exp 



\y\ 



n\g(t)\ 2 *A W)\\ 
yielding a Gaussian function with width given by 



7 



E 



9l 



"1 + (7/2) 2 



-7* _ 2 e -5 t coso; A t + l) , 



which can be express in terms of the spectral density of 
the bath. For the Ohmic case we get 



Markovian channel at zero temperature 

For a Markovian channel at zero temperature coupled 
to a qubit, J — ^J^<j— with er_ the deexcitation operator 
and 7 the decay rate, 0(a, s, s') — J (see Refs. 24, 25 and 
26), and f(a,t) = e" 7 '/ 2 . 

Taking for simplicity a frozen internal dynamics 



H, 



(ch) 

cs — 



0, the equation of motion for the non-unitary 



propagator is 



0_ 
ds 



U (a, s) 



7 

y/jz (a,s)&- - -h- 



U(a, s) 



with initial condition U(a, 0) = I. Here n_ = <7_|_<t_, 
and z(s) = J2\9^ a * e ~ lulxS ^ s a sum over the modes of 



Pit) := Wit)\ 2 = 



2tt _ 
1-e" 7 *. 



00 e" 7 ' - 2e"i*cosa;t + 1 



(7/2) 5 



did 



After collecting our results, the probability distribu- 



tion of entanglement ( 14 ) can be recast into 

_ 2 

e 2 

Pa(x,p) = 



e ^ / e 2 ■ 
5 x 



np(t) 



F{v) 



dy 



with 



F(y) = 5RTr p ch (0) [/ + (e"^ - 1 + \y\ 2 ) n„ + y&+ 
This last integral can be evaluated to finally obtain 
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P G (x,t)=2 



1-P(t) 
v 3 PnP(t) 



x exp 




■ PuP(t) 
Pn 



2jpio| 



Io 



gjpioj / Vl-P(t) _ 1 - pnp(f) 



1/2" 



where lo(z) is the modified Bessel function of zeroth or- 
der. 



Non-Markovian dephasing channel 

We now turn to the general dephasing case in d = d^ ch ^ 
dimensions for which 0(a, s, s') = J = (see refs. 25 

and 26) and consider an arbitrary spectral density. 

Let us first note that in this case 

f(a,t) = e-i T,:J2?li tif« a( - s - s ' )ds ' ds 
is independent of {a( c h)}. Moreover, introducing 

q(t) : = U [ I a(s - s')ds'ds 
Jo Jo 



t 2 f 00 ^ , f smut/2 



I(u>) 



ujt/2 



(15) 
duj. 



we easily get 
with 

iuj\ t 



i 



a\g\- 



The probability distribution for the real variable y(a, t) 



cTgif) is again given by a Gaussian function, 



P v (y,t) 



-\a\ 2 



z 
1 



-Siy-^tcTgit^da 

,2 



■ exp 



y 



with width 



2«(t) := |gWI 2 = E 



WW J' 



-iui x t _ ^ 



The probability distribution of entanglement ( 14 ) for 
this environment simplifies into 



P G (x,t) 



I?(t)Tr J 2 
X 

«!_ 

e 2^xtT 



e -f g(t )TrJ 2 \ 

o I x ^— I dy, 



in which 



If we now restrict to the case in which the coupled 
subsystem is a qubit, i.e., d — 2, the integral can be 
evaluated to obtain the entanglement distribution 



P G (x,t) 



l-p(t) 



xy/l - 4x 2 PnPoo V 2tt| log(l -p(t))\ 



exp 



log' 



2 ( l+Vl-42 



2Iog(l-p(t)) 



exp 



log' 



2 ( l + \/l-4a: 2 piiPon 



2Kpoo 



21og(l-p(t)) 



where we introduced a convenient re-paremeteriztaion of 
timep(i) = l-e- 2TrJ2 9( t ). 



Scaled time for non-Markovian dephasing baths 



The behaviour of the function q(t) for different models 
of the non-Markovian bath is easily found by selecting 
the appropriate spectral density. In the continuum limit, 
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this distribution is assumed to have the general form 

I(w) = Tr- 1 (w/w d ) v e~ w/Wd 

where the cut-off (Debye) frequency uj d is the maximum 
frequency of the bath. 

We may then consider different limiting cases. Adopt- 



ing the standard nomenclature (see Ref. 4), ij -> 0,U(j ^ 
oo defines purely ohmic dissipation, while 77 — > defines 
ohmic dissipation with upper frequency, and 77 — > 1 repre- 
sents superohmic dissipation. Explicit calculation of the 



integral in Eq. ( 15 ) for the different limiting cases leads 
to 



r t 



?(*) = 4 



- tan 1 ir,,l - 1 1 ' 1 — >,J 

7T TtWd 



(purely ohmic) , 
log(l -I- w^i 2 ) (ohmic) , 



TTW d 



log(l + w 2 d t 2 ) 



(superohmic) , 



showing a transition from exponential to algebraic decay 
for the function p(t) = 1 - e -™J 2 iW. 



Convolutionless operator 0(a,s,s') for a 
non-Markovian bath at zero temperature 

In this last section we evaluate the operator 0(a, s, s') 
for a non-Markovian bath at zero temperature, following 
the techniques of Ref. 25. For a CS with frozen internal 
dynamics Hqs — and coupled to the bath trough an 
operator J = ^/7cr_ , indicating amplitude damping, the 
condition for the existence of the operator 

0(a, s, s) = Jit(s, s') 

was shown in Ref. 26 of the main text to be that the 
nonlinear integral equation 

—u(s,s')=^u(s,s') / a(s — s")u(s, s")ds" 

with u(s, s) — 1 admits a solution. 

Remarkably, this nonlinear problem can be exactly 

mapped into a linear one by introducing u(s, s') = 
c(s')/c(s), where 

— c(s) — —7 / a(s — s)c(s')ds . 
os Jo 

If the bath is at zero temperature and has a memory 
kernel of the Ornstein-Uhlcnblcck type 



a(s) =■.■<' 



2 



this last problem takes the form 



d_ 

ds 



which, due to its linearity, can be solved by Laplace trans- 
form methods. After some algebra the result is 



c(s) 



c(0) 



(z_e z+s - z+e z - s ) , 



where 



z± 



are the solutions of the associated characteristic equation, 
and we introduced f2 = (uj/2)y/l — (2jju>dj. Following 
some algebraic manipulations, this result can be cast into 



c(s) = 2c(0)e-^ ; 
from which we recover 
u{s, s') = 



—— sinh (57s) + cosh (Cls) 



<s') = e ^ (s _ s0 m sinh (n S >) + cosh (Sis') 
c(s) ^jf sinh (Sis) + cosh (Sis) 



Since our goal is to construct contribution from this 
particular channel to the exponent of the mean entangle- 
ment x(t), function in Eq. (10) of the main text, we need 
to calculate 

r(s) = / a(s- s')u(s,s')ds', (17) 
Jo 

for this specific choice. Substitution of a(s) and u(s, s') 
gives finally 



r( s ) 



m sinh Qs 



^ sinh Sis + cosh Sis ' 



(18) 



c (s 



-w d (s-s') 



c(s')ds' , (16) thus completing the calculation. 



